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INTRODUCTION 


Analysis  of  radar  detection  performance  in  moving  platforms,  either 
in  a  satellite  or  aircraft,  is  a  continuing  area  of  interest.  The  need 
to  detect  smaller  targets  in  look-down  geometries,  plus  the  expected 
difficulty  of  achieving  deterministic  low  sidelobe  patterns  in  space 
deployable  antennas  and  for  conformal  antennas  on  various  airframes  has 
dictated  a  relook  at  the  moving  target  detection  design  for  such  radar 
systems . 

This  report  describes  a  simulation  analysis  computer  program  which 
examines  the  target  detection  performance  of  space  based  and  airborne 
surveillance  radars  in  look-down  clutter.  The  computer  model  samples 
ground  clutter  with  specified  radar  cross  section  and  grazing  angle 
effect  and  analyzes  clutter  signal  power  intensity  as  an  effective  area 
in  range  and  range-rate  (doppler)  space.  With  included  targets  and 
antenna  parameter  data,  a  computed  target  to  clutter  ratio  plotted  for 
some  range  over  all  range-rates  is  shown  as  a  means  of  measuring  target 
minimum  detectable  velocity.  A  3-D  plot  of  clutter  intensity  as  program 
output  indicates  how  antenna  systems  may  be  compared  with  one  another  in 
terms  of  target  detectability.  The  program  is  capable  of  plotting  in 
3-D  format  the  clutter  intensity  or  clutter  to  target  ratio  over  all 
range  and  doppler  which  is  visible  to  the  antenna.  The  3-D  plots  at  a 
glance  indicate  the  relative  merit  of  given  radar  antenna  systems. 

General  program  parameters  are  developed  with  geometric 
consideration  of  both  a  non-rotating  and  a  more  complex  rotating  earth 
model.  Current  program  limitations  are  discussed  as  well  as  continuing 
and  future  improvement  efforts.  Although  this  program  is  applicable  to 
any  set  of  monostatic  antenna  and  platform  parameters,  a  specific 
example  is  given  as  an  illustrative  example. 


RANGE-DOPPLER  CELLS 

The  problem  of  computing  range-doppler  sample  data  from  a  moving 
platform  is  outlined  in  figures  la  and  lb.  A  range  gated  doppler  radar 
analyzes  the  returned  signal  in  both  time  (range)  and  frequency 
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(doppler).  The  signal  is  processed  by  first  dividing  it  into  a  number 
of  intervals  in  time  (range  gates)  and  then  analyzed  by  passing  it 
through  a  set  of  doppler  frequency  filters.  The  locus  of  points  of 
equal  range  from  a  monostatic  radar  antenna  uhich  intersects  the  earth 
is  called  an  isorange  and  is  a  circle  centered  about  the  platform 
subpoint.  All  such  sampled  ranges  produce  a  series  of  concentric 
circles  on  the  earth  (figure  2a). 


EQUAL  RANGE  CONTOURS  EQUAL  DOPPIER  CONTOURS 

FROM  AN  ELEVATED  PLATFORM  FROM  AS  ELEVATED  MOVING  PLATFORM 


Figure  la.  Figure  lb. 


Unfortunately,  the  locus  of  equal  range-rates  intersecting  the 
earth  is  not  as  easy  to  characterize.  Relative  to  a  moving  platform  the 
locus  of  points  in  space  of  equal  doppler  shift  is  a  cone  whose  vertex 
is  at  the  antenna,  and  whose  axis  coincides  with  the  velocity  vector  of 
the  platform.  Intersection  with  the  spherical  earth  results  in  a 
distorted  conic  section  called  an  isodop.  The  maximum  doppler  component 
of  ground  clutter  is  thus  the  platform  velocity  relative  to  the  earth 
multiplied  by  the  cosine  of  the  minimum  cone  angle  which  intersects  the 
e  ”th  (fip>  e  lb).  The  zero  doppler  component  will  be  at  all  range 
poirtt  wt  :h  lie  orthogonal  to  the  velocity  vector  (neglecting  rotating 


earth  effects) 


The  intersection  of  isoranges  and  isodops  for  a  radar  aystem  with 
contiguous  range  gates  and  doppler  filters  is  sketched  in  figure  2a.  To 
limit  the  processing  required,  the  radar  clutter  returns  may  be  sampled 
at  a  finite  number  of  ranges  equally  spaced  from  the  minimum  to  the 
maximum  expected  range.  The  doppler  data  may  also  be  sampled  by  a 
finite  number  of  doppler  filters  equally  spaced  between  maximum  and  zero 
doppler.  The  minimum  width  of  the  range  gates  may  be  taken  down  to  the 
range  resolution  limit  of  the  radar  system.  Similarly,  the  doppler 
filter  bandwidth  minimum  limit  depends  on  the  radar  system 
characteristics.  The  intersection  of  N  lsoranges  with  M  isodops  and 
their  respective  resolution  widths  define  NxM  range-doppler  cells  (see 
figure  2b).  It  is  from  these  cells  that  the  simulation  clutter  data  is 
sampled  and  computed. 


Figure  2a.  Equal  Range  and  Eaual 
Doppler  Map  on  Earth's  Surface 


Figure  2b.  Sampled  Resolution  Limited  Map 


GEOMETRIC  ANALYSIS 

Figure  3  shows  a  moving  platform  radar  antenna  at  a  distance  H 
above  the  earth's  surface  In  a  right  hand  coordinate  system.  A 
geometric  analysis  of  figure  3  yields  the  equations  necessary  to 
calculate  the  range~doppler  cells  and  their  effective  clutter 
contribution  to  the  radar  echo. 

Any  point  in  the  far  field  of  an  antenna  can  be  represented  by  a 
pair  of  direction  cosines  Tx  and  ,  where  the  relationship  with 
spherical  coordinate  angles  is: 

-  C08(Ax)  -  sin(e)cos(<t> ) 

T^.  -  co8(Ay)  -  sin(6)sin(<J>) 

It  is  convenient  to  place  the  moving  antenna  platform  at  the  center 
of  the  coordinate  system  with  the  velocity  vector  V  coinciding  with  the 
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Y  axis  and  the  Z  axis  pointing  up.  The  reference  to  a  point  on  the 
ground  may  be  given  in  terns  of  a  vertical  elevation  angle  C  and  an 
azimuthal  scan  angle  S  measured  from  7.  The  direction  cosines  are 
simply  related  by: 

Tx  -  sin(C)sln(S) 

Ty  -  sln(C)cos(S) 

From  figure  4a  the  elevation  angle  of  the  horizon,  which  is  the 
complement  of  the  minimum  horizon  cone  angle  and  the  range  to  the 

visible  horizon  R^y  are: 

sln((1nax>  *  coa<Vnin>  *  R*/(RefH) 

Umax  -  ( (RefH) 2  -  Re2)% 

T 

B 

i  * 

R« 

Jl. 


Figure  Ua. 

Geometry  at  the  Horizon 


Figure  Ub. 

Geometry  at  Intermediate  Ranges 


Similarly,  for  any  given  range  R(I)  corresponding  Co  some  range  bln  1, 
figure  4b  and  the  law  of  cosines  results  In  the  elevation  angle  C(I): 


C(I)  ■  arccos 


R2(I)  +  (R  +H)2  -  R2 
e  e 


2R(l) (R  +H) 
e 


The  grazing  angle  G( I)  at  a  particular  range  R(I)  Is  given  by 
figure  4b  and  the  law  of  sines: 


G(I)  ■  arccos 


sin(C(l)) 

R 


Consider  now  the  equally  spaced  sampled  range  and  doppler  bins  of 
figure  2b.  The  distance  between  the  N  range  bins  Is  determined  by 
dividing  the  total  range  difference  by  the  number  of  range  gates  less 
one: 


Dr  ■  (Rmsur*  H)/(N-1) 

The  range  of  the  1th  (1-1  to  N)  bln  Is  calculated  from: 

R(I)  -  (I-l)Dr  +  H 

where  the  nadir  Is  at  R(1)«H  and  the  horizon  Is  at  R(I~N)-Rmax. 

To  avoid  the  problems  of  clutter  cell  area  on  the  horizon,  for  the 
purposes  of  the  simulation  the  range  difference  Is  divided  by  N  Instead 
of  N-l .  This  has  the  effect  of  decreasing  Dr  slightly  and  placing  the 
last  range  bin  R(N)  a  small  distance  back  from  the  horizon.  If  the 
range  resolution  of  the  radar  is  Ar,  from  figure  4b  the  effective  ground 
range  coverage  for  small  Ar  Is  Ar/cos(G(I)).  Of  course  N  must  be  chosen 

so  that  D  Is  significantly  larger  than  this  value  to  prevent  range 

r 

overlapping  In  the  range  bins  of  interest. 

The  M  equally  spaced  doppler  bins  are  equivalently  computed  from 
the  radial  velocity  components,  since  the  two  way  radar  doppler  shift  is 
proportional  to  the  radial  velocity  seen  at  the  antenna  platform  and  the 
wavelength  of  the  radar  operating  frequency: 
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F  -  2V  /X 
d  r 


where  Che  radial  velocity  component  is  related  to  the  platform  velocity 
(for  a  stationary  earth): 

V  -  V  cos (A  ) 
r  p  y 

Then  the  frequency  difference  between  the  M  equally  spaced  doppler 
filters,  and  the  cone  angles  which  yield  the  M  isodops  are: 

D  *  2V  cos(A  )/mx 
f  p  ymin 

Ay(J)  «  arccos(JDf  X/2Vp)  J"1  to  M 

If  the  doppler  filter  width  is  Af,  then  the  angles  which  yield  the  N 
Isodops  of  width  Af  are  between  Ay(J)  and: 

A  (J)  ■  arccos(( JD,+Af )X/2V  ) 

yd  f  p 

A  critical  factor  arises  at  this  point  in  the  calculations.  Since 
a  cosine  function  has  its  most  rapid  change  at  large  angles,  most  of  the 
equally  spaced  doppler  bins  will  be  crowded  at  cone  angles  greater  than 
45  degrees  away  from  the  velocity  vector.  By  the  same  reasoning,  the 
ground  coverage  due  to  these  angles  is  larger  for  smaller  cone  angles 
(i.e.  close  to  the  plane  of  motion  near  the  horizon).  In  effect  the 
area  of  the  range-doppler  bins  of  figure  2b  becomes  larger  as  the  plane 
of  motion  is  approached. 

Now  the  clutter  intensity  in  the  radar  echo  will  be  a  function  of 
the  area  of  the  range-doppler  cells  in  figure  2b.  From  here  on  a 
simplifying  assumption  must  be  made  concerning  the  shape  of  the 
range-doppler  cells.  As  long  as  the  resolution  width  of  the  range  bins 
and  doppler  bins  is  small  compared  to  the  distance  between  them,  each 
cell  can  be  approximated  as  a  rectangular  area.  The  radial  component  is 
a  simple  function  dependent  on  the  grazing  angle  and  range  resolution, 
but  the  lateral  component  requires  a  somewhat  more  complicated  geometric 
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analysis  of  the  vectors  which  define  the  isodops  and  their  equivalent 
resolution  widths  at  each  range. 


CLUTTER  INTENSITY  DIAGRAM 

Consider  the  case  of  continuous  scatterers  uniformly  distributed  on 
the  earth's  surface  and  a  radar  with  infinitely  fine  resolution.  If  the 
radar  is  unambiguous  in  both  range  and  doppler,  the  clutter  scatter 
signal  strength  C(r,r)  will  be  a  function  of  at  least  four  components. 
First  is  the  1/R  two  way  path  loss  attenuation  factor  from  the  radar 
range  equation.  Second  is  the  radar  cross  section  aQ  ,  of  the  ground 
scatterers.  Several  models  are  available  which  describe  ground  RCS  as  a 
function  of  terrain  type  and  grazing  angle.  The  third  component  is  the 
two  way  relative  antenna  gain  as  a  function  of  directional  cosines  Tx 
and  .  The  directional  cosines  correspond  to  those  angles  which  yield 
r  and  r  at  the  earth's  surface.  Finally,  the  area  of  the  range-doppler 
cells  must  be  computed  which  will  determine  the  number  of  discrete 
scattering  centers,  providing  the  integrated  clutter  signal  power. 

Presuming  for  the  moment  an  omnidirectional  antenna  with  unity  gain 
and  uniformly  distributed  range-doppler  scatterers,  the  clutter 
intensity  diagram  takes  on  the  smoothly  shaded  appearance  of  figure  5a. 
Range  R  is  plotted  along  the  ordinate  axis  and  the  direction  cosine 
angle  Ay  is  plotted  on  the  abscissa.  Dark  shading  represents  regions  of 
high  clutter  echo  intensity  (integrated  radar  return  signal  power). 

The  distinctive  "U"  shape  of  the  expected  clutter  return  of  figure 
5a  is  strictly  a  function  of  platform  altitude  and  earth  radius. 
Outside  this  shaded  area  is  a  clutter  free  region  where  no  combination 
of  range  and  cone  angle  intersect  the  earth's  surface  and  in  this  region 
a  target  would  not  compete  with  surface  clutter.  Directly  beneath  the 
antenna  platform  corresponds  to  a  zero  degree  elevation  angle  and  90 
degree  cone  angle,  with  range  equal  to  platform  altitude  (figure  lb). 
Progressing  outward  in  range  one  reaches  the  radio  horizon,  which  is  the 
range  of  the  visible  horizon  from  the  elevated  antenna  platform. 
Greater  ranges  would  correspond  to  range-doppler  cells  which  do  not 
intersect  the  earth.  Range  excursion  along  the  central  vertical  of 
figure  5a  at  a  constant  90  degree  cone  angle  is  a  vertical  scan  in 


elevation  orthogonal  to  the  plane  of  aotion,  while  travel  along  the  ”0" 
edge  la  a  vertical  scan  in  the  plane  of  aotion.  A  depiction  of  the 
relationship  between  a  pencil  beaa  antenna  footprint  and  the  clutter 
Intensity  dlagraa  is  shown  in  figure  Sb.  Notice  that  near  the  edge, 
close  to  the  plane  of  aotion,  the  intensity  rises  due  to  the  Increased 
area  size  of  the  range~doppler  cells.  Since  the  cell  area  is  larger, 
effectively  aore  point  scatterers  are  covered  per  cell  resulting  in  a 
greater  power  intensity. 
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Figure  5a»  Figure  Sb. 

Clutter  Intensity  Diagram  C(r,r)  Relationship  Between  a  Pencil  Beam  Antenna 

Footprint  and  the  Clutter  Intensity  Diagram 
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Mow  consider  a  somewhat  more  realistic  antenna  with  finite  range 
resolution  &r  and  finite  doppler  resolution  A r  corresponding  to  discrete 
noncontiguous  range-doppler  cell  scatterers.  The  clutter  scatter 
function  for  each  cell  is  depicted  in  figure  6.  This  discrete 
characterization  of  the  distributed  returns  C(Ar,Ar)  for  the  ^ 
specific  range  bin  and  the  specific  doppler  bin  is  used  in  the 

computer  program  simulation.  The  clutter  scatter  function  of  each  cell 
is  an  integral  over  the  range  and  doppler  cell  widths. 


C.  (Ar,Ar) 
i » J 


r  +Ar/2  r  +Ar/2 

J  .J 

r.-Ar/2  r.-Ar/2 

i  J 


C(r.,fJ)dri4rj 


Figure  6. 

Clutter  Intensity  Diagram 
Discrete  Scatterer  Case  C(Ar,Ar) 


For  an  unambiguous  range  gated  pulse  doppler  radar  with  a  known  far 

field  antenna  pattern  weighting  (including  sidelobes)  as  a  function  of 

directional  cosines,  one  only  needs  to  multiply  the  discrete  clutter  by 

the  weighting  function  ANT(T  ,T  )  to  obtain  the  clutter  intensity 

x,  y 

function  of  the  antenna. 
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DIRECTING  THE  ANTENNA 

The  qu Melon  of  directing  the  antenna  mein  beam  toward  some 
specified  point  on  the  earth's  surface  must  be  addressed.  In  a  phased 
array  this  is  accomplished  by  electronic  beam  steering,  which  simply 
redefines  the  weighting  function  ANTC^.T^).  Alternatively,  the  antenna 
platform  may  be  mechanically  steered  in  roll,  pitch,  and  yaw  as  defined 
in  the  right  hand  coordinate  system  of  figure  7. 

2 


Figure  7. 

Coordinate  Transformation  Axes 


Mechanical  steering  is  accomplished  by  applying  the  Jacobian 
matrix: 


cos(pit)cos(yav) 


-cos(pit)sin(yav) 


T 


x 


rn  t 

A 

y 


sin(rol)sin(pit)cos(yaw) 

+cos(rol)sin(yav) 


-sin(rol)sin(pit)sin(yav) 
+cos ( rol ) cos ( yaw) 


T 

y 


(From  Mathematical  Handbook  for  Scientists  and  engineers,  Korn  and  Korn 
2nd  ed . ,  McGraw  Hill,  1968) 


where  the  primed  coordinates  are  in  tha  staared  fraaa  of  reference. 
These  new  directional  cosines  are  input  to  the  aodifiad  weighting 

function  ANT(T  ',T  ').  The  antenna  has  been  effectively  steered  toward 

*  y 

the  specified  patch  of  earth. 


ROTATING  EARTH  MODEL 

Up  to  this  point  the  effects  of  the  rotating  earth  have  been 
neglected.  This  assumption  is  valid  for  airborne  antenna  platforms 
which  are  in  a  rotating  earth  frame  of  reference,  but  the  motion  of 
space  based  radar  platforms  is  Independent  of  earth  rotation  and  the 
earth  will  provide  a  component  of  doppler  velocity.  Consider  a  platform 
stationary  in  space  with  the  earth  rotating  beneath  as  sketched  in 
figure  8. 


i 


Figure  8. 

Rotating  Earth  Doppler  Components 


Clearly  for  an  orbiting  platform,  the  earth's  doppler  velocity 
component  will  be  negative  when  looking  in  the  eastern  hemisphere, 
positive  in  the  western  hemisphere,  and  zero  along  the  central  meridian. 
Of  course,  to  the  earth  component  radial  velocity  must  be  added  the 
radial  component  due  to  the  orbital  velocity  of  the  satellite.  Note 
Chat  the  line  of  zero  total  doppler  velocity  will  no  longer  be  precisely 
orthogonal  to  the  velocity  vector  (except  in  an  equatorial  orbit).  In 


13 


general  for  a  moving  platform  In  circular  orbit,  the  radial  doppler 
velocity  component  depends  on  altitude,  orbit  Inclination,  and  the 
grazing  angle  of  the  ground  clutter  patch.  In  the  simulation  of  a 
range-doppler  cell  radar  system  the  total  radial  velocity  component  from 
each  of  the  individual  cells  is  given  by: 

Vr(I,J)  *  VpCos(Ay)  -  Vecos(G)cos(lat)cos(scan) 

scan  -  arccos(cos(Ay)/8ln(C))  -  arccos(cos(inc)/cos(lat)) 

where  A  Is  a  function  of  I  and  J;  G  and  C  are  functions  of  I;  "inc"  is 

y 

orbital  inclination;  "lat"  Is  instantaneous  satellite  sub-latitude;  and 
Ve  is  the  equatorial  rotational  velocity  of  the  earth.  The  term  "scan" 
is  essentially  the  azimuthal  scan  angle  measured  from  the  satellite's 
instantaneous  velocity  vector. 

The  maximum  total  doppler  velocity  component  will  be  at  the  horizon 
In  the  plane  of  motion.  (Actually,  depending  on  the  position  of  the 
velocity  vector,  the  rotating  earth  may  change  this  maximum  component 
direction,  but  for  circular  orbits  up  to  geosynchronous  altitudes  the 
effect  is  small.) 

^rmax  "  ^pcos^ymin^  “  Vecos(inc) 

The  velocity  differential  between  each  of  the  M  equally  spaced  doppler 
bins  is  thus: 


V  ■  V  /M 
diff  rmax  ' 


and  the  doppler  velocity  of  the  J  doppler  bin  is 


JVoax  ■  V’J) 


Solving  for  the  cone  angles  Ay(I,J)  which  yield  the  desired  equally 
spaced  doppler  bins  results  in  a  considerably  more  complex  function  than 
for  the  nonrotating  earth  model  because  of  the  combined  I,J  dependence. 
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where: 


0V1  •  (7  cos(A  )  -  V  cos(inc))J/M 

p  vtnin  e 

DV2  •  Vp  -  VgCoaUnO/cosU^j,  ) 

0V3  -  (cos2(lat)  -  coa2(inc))l5Ve/co8(Ayiniri  ) 

SQV  »  [dV22  +  DV32  -  DVl2/sin2(C(I))]  ^ 

If  Che  doppler  filter  bandwidth  is  Af ,  then  the  angles  which  yield 
the  N  lsodops  of  width  Af  are  between  A^I,J)  and: 

AydfI,J^  "  Av^1, J^DV1=DV1  +  AfX/2 

Care  oust  be  exercised  when  applying  these  equations  to  take  into 
account  the  direction  of  orbit  and  whether  the  antenna  is  pointed  toward 
the  northern  or  southern  hemisphere.  Por  the  nonrotacing  earth  case 
(Ve-0),  or  for  the  satellite  subpoint  over  a  pole  (inc-90,  lat-90)  the 
equations  reduce  to  the  simpler  form  previously  derived. 


THE  CLUTTER  SIMULATION  PROGRAM 

The  program  CLUTR  is  designed  for  either  a  space  based  radar 
antenna  in  circular  orbit  or  for  an  airborne  platform  in  level  flight. 
The  antenna  far  field  weighting  input  data  resides  in  a  file  which 
either  has  been  previously  generated  from  PAAS  (Parametric  Antenna 
Analysis  Software)  or  from  an  actual  antenna  data  cape. 

The  Inputs  for  a  circular  orbiting  space  based  moving  platform 
radar  are: 

(1)  The  altitude  above  the  earth's  surface. 

(2)  The  range  and  doppler  bin  resolution  size. 

(3)  The  ra<*nr  operating  frequency. 

(4)  The  orbit  Inclination  and  platform  sub-latitude. 

(3)  The  direction  of  orbit,  east  or  west. 

(6)  The  main  beam  pointing  quadrant,  north  or  south. 


(7)  Th«  antenna  roll,  pitch,  cod  yaw  pointing  angles. 

(8)  The  target  radar  cross  section. 

(9)  The  far  field  antenna  weighting  values  ANT(l^,7y). 

Inputs  for  an  airborne  platform  are  the  ease  except  Items  (4),  (S),  and 
(6)  are  omitted  and  an  input  ground  speed  Is  required.  Figure  9  is  a 
block  dlagraa  of  the  basic  components  of  the  program. 


Figure  9. 

Program  CLUTR  Block  Diagram 

The  data  from  two  auxiliary  programs  are  required  for  Inputs  to  the 
main  program  CLUTR.  This  is  to  insure  that  the  radar  main  beam  falls 
within  one  of  the  sampled  range-doppler  cells.  Program  RANGE  provides  a 
listing  of  range  bin  numbers  with  range  R(I)»  elevation  angle  C(I),  and 
grazing  angle  G(I)  for  each  range  bin  I.  Given  a  desired  grazing  angle, 
for  example,  the  range  bin  number  yielding  the  nearest  grazing  angle  is 
selected.  This  grazing  angle  and  corresponding  elevation  angle  value 
must  be  entered  iuco  CLUTR  and  will  center  the  main  beam  in  range  I. 
Program  SCAN  provides  a  listing  of  azimuthal  scan  angle  S(I,J)  and  cone 
angle  Ay(I,J)  for  each  doppler  bin  J  at  the  range  I.  The  value  S(I,J) 
centers  the  main  beam  in  doppler  bin  J. 
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EXAMPLE  CLUTTER  PROGRAM  RUN 


The  utility  of  the  simulation  is  best  demonstrated  with  an  example 
case.  Consider  first  an  omnidirectional  antenna  (i.e.  all  and  T 

antenna  weighting  values  are  unity)  with  the  antenna  platform  in  a  4000 
nautical  mile  high  circular  orbit.  Samples  of  128  range  bins  and  128 
doppler  bins  are  taken,  resulting  in  128  range-doppler  cells.  The  3-D 
sampled  clutter  intensity  diagram  output  is  shown  in  figure  10.  Notice 
that  figure  10  is  just  one  half  of  figure  6,  with  the  inter-cell  spaces 
compressed.  The  ragged  edge  is  an  artifact  of  the  sampling  process,  as 
no  samples  are  taken  exactly  in  the  plane  of  motion,  but  the  distinctive 
"U"  shape  is  clearly  evident. 


Figure  10. 

3-D  Omnidirectional  Antenna  Clutter  Inters itv 


Now  place  Che  antenna  aboard  an  orbiting  platform  with  the 
following  observation  characteristics: 

Altitude  -  4000  nautical  miles 

Inclination  -  60  degrees 

Sub-latitude  Point  -  20  degrees 

Main  Beam  Scan  -  True  north  (000  degrees  true) 

Desired  Main  Beam  Grazing  Angles  -  5  degrees  and  30  degrees 
Target  RCS  -  20dBsm 


From  the  auxiliary  program  output  RANGE,  the  grazing  angle  closest 
to  30°  is  30.04^  ,  corresponding  to  an  elevation  angle  of  23 • 5 sP  at 
range  bin  number  55.  Similarly  from  the  same  output  listing,  a  grazing 
angle  of  4.84°  is  in  range  bin  number  115  with  an  elevation  of  27. 4^  . 
The  antenna  is  steered  by  providing  the  elevation  angle  as  a  roll  or 
pitch  angle. 
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Sample  Output  Propram  RANGE 
5°  Grazing  Angle  Case 


At  this  point  in  the  orbit,  the  velocity  vector  heading  relative 
to  true  north  is: 

HDG  ■  arc8in(cos(inc)/co8(lat)) 

-  32.15  degrees 
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Thus  Co  scan  to  a  true  heading  of  000  degrees  a  steering  yaw  of  32.1S0 
is  required.  From  the  auxiliary  program  output  SCAN,  the  nearest  scan 
angle  for  the  55  th  range  bin  is  32.718°  in  the  103rc^  doppler  bin;  for 
the  115th  range  bin  it  is  32.232°  in  the  119th  doppler  bin. 
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Sample  Output  Program  SCAN 
5°  Grazing  Angle  Case 


Figure  12  is  a  3-D  plot  of  the  antenna  pattern  weighted  clutter 
intensity  diagram  for  Che  5°  grazing  angle  case.  The  main  beam  clutter 
region  is  clearly  defined  (centered  at  range  bin  number  115  and  doppler 
bin  number  119),  although  it  is  smeared  in  range  near  the  horizon.  At 
low  grazing  angles  such  as  this,  even  if  the  main  beam  has  relatively 
narrow  angular  beamwldth,  It  will  cover  many  range  clutter  cells  due  to 
the  curvature  of  the  earth.  Again  the  "U"  shape  similar  to  figures  5a 
and  10  is  apparent,  with  most  of  the  clutter  along  the  doppler  edge. 
Even  though  the  main  beam  is  pointed  well  out  in  range  and  the  sidelobes 

L4 

are  very  low,  there  is  some  close  in  clutter  do  to  (1)  the  1/R  factor 
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which  has  much  less  effect  for  small  R,  and  (2)  the  grazing  angle  effect 
which  increases  the  ground  patch  clutter  RCS  for  near  ranges. 


A  similar  plot  in  figure  13  is  for  the  30°  grazing  angle  case. 
Here  the  beamwldth  on  the  ground  is  relatively  narrow  and  the  sidelobes 
still  contribute  to  the  near  range  clutter,  although  performance  is 
better  because  there  is  not  such  a  large  difference  between  the  range  of 
the  main  beam  clutter  and  the  nadir. 


ANTENNA  WITH  BLOCKAGE 

Consider  the  same  PAAS  generated  antenna  pattern,  only  this  time 
with  a  5%  central  blockage.  The  far  field  pattern  is  shown  in  figure 
14.  Notice  that  the  blockage  has  significantly  degraded  the  aldelobe 
structure  so  that  this  antenna  would  be  expected  to  have  degraded  off 
axis  clutter  suppression. 


Figure  14. 

Uo  dB  Taylor  Far  Field  Antenna  Pattern 
With  5%  Central  Blockage 
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Figure  15  is  a  3-D  plot  of  the  clutter  intensity  diagram  for  the  5 
grazing  angle  case  with  the  blockage  antenna.  The  general  form  is,  as 
expected,  the  same  as  figure  12  but  a  considerable  amount  of  sidelobe 
clutter  has  been  introduced.  Figure  16  is  similarly  a  degraded  version 
of  figure  13. 


iv'i'H  4:fo3Ti,v 


Figure  IS. 

3-D  Clutter  Intensity  Plot  With  Blockage 
S°  Grazing  Angle 
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Figure  1 6. 


30  Grazing  Angle 


MINIMUM  DETECTABLE  VELOCITY 


The  target  to  clutter  ratio  output  of  CLUTS.  for  the  55th  range  bin 
over  all  doppler  bins  is  plotted  in  figure  17a  and  expanded  in  figure 
17b.  Output  for  the  115th  range  bln  is  plotted  in  figures  18a  and  18b. 


NMM(V«  »»M  MMIVB 


Figure  17a. 

Target  to  Clutter  Ratio  Versus  Doppler  Bin 
30°  Grazing  Angle 
Range  Bin  55 
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Figure  17b. 

Target  to  Clutter  Ratio  Versus  Doppler  Bin 
30  Grazing  Angle 
Range  Bin  55 


-  25  - 


MMira  iin  Nun»u 


Figure  l8a. 


Target  to  Clutter  Ratio  Versus  Doppler  Bin 
5°  Grazing  Angle 
Range  Bin  115 
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Notice  that  the  minimum  plot  value  falls  in  the  103rd  doppler  bin 
for  the  55^  range  bin  of  figure  17b.  This  doppler  bin  corresponds  to 
the  radar  return  of  a  target  which  is  stationary  with  respect  to  the 
earth.  The  rise  to  an  infinite  target  to  clutter  ratio  at  doppler  bin 
112  corresponds  to  cone  angle  (target  doppler)  and  elevation  angle 
combinations  which  do  not  intersect  the  earth  and  hence  the  target  does 
not  compete  with  any  clutter  (clutter  free  region).  Similarly,  in 
figure  18b  the  stationary  target  minimum  value  falls  in  doppler  bin  119. 

The  fixed  clutter  echoes  with  which  the  desired  target  signal  mu6t 
compete  are  those  included  within  the  same  radar  resolution  cell  as  the 
target,  or  those  which  enter  the  receiver  via  the  antenna  sidelobes.  To 
interpret  the  above  plots  in  terms  of  minimum  detectable  velocity,  note 
that  the  minimum  value  of  the  target  to  clutter  ratio  corresponds  to  a 
stationary  target  (relative  to  the  earth)  located  in  the  same 
range-doppler  cell  as  the  clutter  patch  that  the  main  beam  is  pointed 
to.  This  cell  location  is  set  to  zero  relative  velocity  between  ground 
and  target.  Thus,  in  order  to  detect  a  stationary  target  a  hypothetical 
receiver  would  require  a  detectability  of  better  than  -20dB  at  a  30° 
grazing  angle  and  better  than  -12dB  at  a  5°  grazing  angle.  However  as 
the  target's  radial  velocity  increases,  detectability  performance 
increases  markedly.  Each  doppler  bin  number  is  equivalent  to  a  radial 
velocity  difference  of  34.3  knots  (another  output  of  SCAN),  which  is 
constant  over  all  range  bins.  At  a  30°  grazing  angle,  with  a  target 
velocity  of  +64.6  knots  (two  doppler  bins)  the  receiver  detectability 
requirement  jumps  to  a  more  reasonable  value  of  25dB.  Similarly  at  a  5° 
grazing  angle,  a  target  velocity  of  +64.6  knots  increases  the 
requirement  to  32dB. 

A  composite  of  all  128  target  to  clutter  ratio  range  bins  produces 
the  3-D  clutter  to  target  ratios  of  figures  19  and  20.  As  clutter  is 
the  dominant-  factor  and  the  target  RCS  is  non  fluctuating,  the  3-D 
clutter  to  target  ratio  plots  are  near  copies  of  their  respective  3-D 
clutter  plots. 
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LIMITATIONS 

The  waveform  characteristics  are  not  Included  in  the  simulation 
model.  There  are  two  consequences  of  this  limitation.  First,  the  total 
energy  transmitted  in  the  radar  range  equation  is  unknown  and  as  a 
result  there  can  be  no  absolute  level  set  on  probability  of  detection  of 
any  given  radar  system;  only  the  ratio  of  target  to  clutter  return 
signal.  Second,  the  particular  shape  of  the  waveform  determines  the 
resolution  size  (above  that  already  determined  by  the  physical  antenna 
parameters)  and  the  ambiguity  function.  The  true  radar  return  signal 
will  be  the  convolution  of  the  clutter  intensity  diagram  and  the 
ambiguity  function. 

The  minimum  detectable  velocity  is  computed  radial  to  the  antenna 
platform,  and  as  such  the  actual  velocity  of  the  target  over  the  ground 
may  be  much  larger  in  specific  cases.  This  is  especially  true  for  large 
grazing  angles  where  the  target  is  illuminated  from  overhead.  Dividing 
by  the  cosine  of  the  grazing  angle  will  in  some  sense  result  in  a 
minimum  ground  speed  in  the  radial  plane. 

The  simulation  presumes  a  "snapshot"  radar  image,  and  then  only  at 
designated  sampled  range-doppler  cell  positions.  The  intersection  of 
the  beam  with  the  ground  for  a  moving  platform  is  a  dynamic  process 
(even  more  so  for  a  moving  target)  which  in  a  real  system  must  be 
considered  over  an  interval  of  time.  Additionally,  the  target  RCS  is 
assumed  nonfluctuating  and  independent  of  aspect  angle. 

The  simulation  program  orbit  is  presumed  circular.  Many 
surveillance  space  radar  applications,  however,  are  in  highly  elliptical 
orbits.  A  minor  program  modification  will  handle  elliptical  orbits  if 
the  orbit's  radial  and  tangential  velocities  are  known. 

To  completely  describe  the  far  field  antenna  pattern  of  a  radar 
system  requires  an  enormous  amount  of  computer  memory  space.  The 
pattern  can  only  be  sampled  at  a  finite  number  of  discrete  points  in  T 
space.  The  quality  of  the  program  results  depends  to  a  large  extent  on 
the  fidelity  of  the  far  field  pattern  and  the  extent  to  which  the 
sldelobe  structure  is  detailed. 
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FUTURE  DEVELOPMENT  EFFORTS 


The  program  CLUTR  is  in  a  continuous  stste  of  change  and 
improvement.  Current  goals  lead  to  inclusion  of  the  radar  waveform  and 
eventual  convolution  with  the  ambiguity  function.  This  will  result  in  a 
true  target  to  clutter  ratio  prediction  for  a  given  radar  platform  and 
will  allow  quantitative  analysis  of  the  system.  The  present  program  is 
only  applicable  to  monostatic  radars  with  no  ambiguities,  but  the 
increasing  importance  and  development  of  bistatic  radar  systems  will 
dictate  program  expansion  to  bistatics  capabilities  as  well.  Ongoing 
future  developments  seek  to  expand  the  utility  of  the  simulation  with 
greater  accuracy  and  a  wider  range  of  applications. 


CONCLUSION 

A  major  problem  in  designing  a  radar  system  with  look-down 
capability  is  determining  target  detectability  in  the  presence  of  ground 
clutter.  The  simulation  analysis  program  is  useful  for  the  designer  of 
a  range  gated  pulse  doppler  to  determine  the  clutter  levels  with  which  a 
target  in  a  given  cell  must  compete  for  detection.  The  results,  while 
only  a  model  and  subject  to  certain  limitations,  should  be  most  useful 
for  comparing  various  candidate  radar  systems  in  the  early  concept 
stages.  Considerable  engineering  savings  in  time  and  cost  may  be 
realized  by  the  early  elimination  of  inferior  candidate  antenna/platform 
system  designs.  The  procedure  described  in  this  report  should  lead  to 
more  efficient  radar  system  tradeoff  analyses  for  the  design  engineer. 
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Listing  of  Program  CLUTR 
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Listing  of  Program  SCAN 
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APPENDIX  -  IV 


Glossary 


ANT(^,Ty)  Antenna  pattern  weighting  function. 

Ax  Polar  angle  from  X  axis. 

Ay  Polar  angle  from  Y  axis. 

Ayj  Polar  angle  from  Y  axis  with  doppler  increment  included. 

Aymin  Minimum  Y  axis  angle  which  intersects  the  horizon. 

C(I)  Elevation  angle  as  a  function  of  range. 

cmax  Maximum  elevation  angle  which  Intersects  the  Earth. 

C(r,r)  Clutter  intensity  scatter  function,  dependent  on  range  and  doppler. 
Df  Frequency  separation  of  consecutive  doppler  bins. 

Dr  Range  separation  of  consecutive  range  bins. 

Fd  Two  way  doppler  shift  of  radar  echo. 

G(I)  Grazing  angle  as  a  function  of  range. 

H  Height  of  radar  platform  above  the  Earth's  surface, 

inc  Satellite  orbit  inclination  angle. 

Isodop  Locus  of  equal  doppler  points  on  the  Earth’s  surface. 

Isorange  Locus  of  equal  range  points  on  the  Earth's  surface, 
lat  Instantaneous  satellite  subpoint  latitude. 

R(I)  Range  as  a  function  of  range  bln  number. 

Re  Radius  of  the  Earth. 

R max  Range  of  the  most  distant  range  bin. 

Tx  Directional  cosine  from  X  axis. 

Ty  Directional  cosine  from  Y  axis. 

Vdiff  Ve1oc1ty  differential  between  successive  doppler  bins 
V@  Equatorial  velocity  of  the  rotating  Earth. 

V  r  Radial  velocity  component  of  a  target. 

V___w  Maximum  relative  radial  velocity  of  the  Earth's  surface. 

rmax 

Vp  Radar  platform  alr/orbltal  velocity. 

Af  Doppler  filter  bandwidth. 

Ar  Range  cell  resolution  width. 

A4-1 


MISSION 
of 

Rome  Air  Development  Center 

HAVC  plans  and  executes  research,  development,  tedt  and 
detected  acquisition  urograms  in  support  of  Command,  Control 
Communicatio ns  and  I tielligence  (C3 1)  activitie d.  Technical 
and  engineering  support  within  area*  of  technical  competence 
ii  provided  to  ESP  Program  Offices  (PCs)  and  other  ESP 
elements.  The  principal  technical  million  arejas  aae 
communication d,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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